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For stationary axially symmetric spacetimes we find a simple expression for the Lense-Thirring 
precession in terms of the Ernst potential. This expression is used to compute, in the weak field 
approximation, the major non-spherical contributions to the precession of a gyroscope orbiting 
the Earth. We reproduce previously known results and give a new estimation for non-spherical 
contributions. 
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I. INTRODUCTION 

The Gravity Probe B satellite launched in 2004 is a 
test of the general theory of relativity sensitive to off- 
diagonal components of the metric tensor resulting from 
the Earth rotation. The satellite contains a set of gy- 
roscopes [1[ that, according to Einstein theory of grav- 
ity, over the course of a year will precess about 6.606 
arcsec/year(0.0018 degrees/year) in the orbital plane, 
an effect known as geodetic precession and 39 milliarc- 
sec/year(0. 000011 degrees/year) in the plane of the Earth 
equator, an effect known as Lense-Thirring precession 
00. The experiment is now in its final phase of data 
analyze [3| and the results are intended to be published 
in the present year The Lense-Thirring part of the 
precession is the one related directly to the off diagonal 
components of the metric tensor and so its experimental 
verification will test the Einstein theory for gravitation. 

The aim of this paper is to relate the expression for 
the Lense-Thirring precession with the mathematically 
sound theory of multipolar moments for axially symmet- 
ric spacetimes due to Geroch [j| and Hansen @ and to 
stimate the multipolar corrections based on some earth 
models. In the weak field approximation we express our 
formulas in terms of Thorne moments Q @| using an har- 
monic coordinate system. This is useful when we consider 
spacetimes whose interior solutions are not known. The 
generalization to non linear fields does not seem difficult 
and can be made using the full Thorne metric @ ■ 

The paper is outlined as follows, in Section [Til we con- 
sider an axially symmetric stationary spacetime and use 
an appropriate tetrad basis to obtain an exact expression 
for the Lense-Thirring precession. By using the Thorne 
metric we specialize this expression for the case of a weak 
gravitational field. In Section IIIII we discuss the pre- 
cession for the Gravity Probe B experiment using three 
known models for the Earth metric and a model based 



in the exact solution of Einstein equation for a rotating 
mass with a quadrupole moment. In particular, we dis- 
cuss multipolar corrections. In Section IV we summarize 
and discuss our results. 



II. THE VECTOR tl LT 

Along the paper we use units such that c = G = 1 and 
metric signature (+,+,+,—). The symbol <> denotes 
the symmetric and trace free part of a tensor. Ai is a 
shorthand notation for a ai ...a ai . Greek letters run from 
1 to 4 and Latin from 1 to 3. 

From the weak field approximation of General Relativ- 
ity and the fact that the spin of an orbiting gyroscope is 
Fermi- Walker transported along its worldline [lCj] we can 
show [2] [3] that the gyroscope spin will precess with an 
angular velocity f2 that can be separated in three parts 
according to their physical origin, 



where 
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The expression ([T]) is general and it is valid for any gy- 
roscope describing a timelike world line. The first term 
is known as the Thomas precession, it depends on the 
gyroscope three-acceleration, a, as well as, on the three- 
velocity, v. It is null when the gyroscope moves along a 
geodesic. The second term is known as the geodetic or de 
Sitter precession and is related to the Newtonian poten- 
tial U. The third term is the angular velocity associated 
to the Lense-Thirring effect. It is related to non-diagonal 
part of the metric, h = g^jej, where denotes the spa- 
tial part of a tetrad base. An important characteristic 
of this term is that it does not depend on the gyroscope 



2 



velocity along the wordline. This fact suggests to study 
the gyroscope precession when it is at rest relative to a 
far observer. Since v is null in such a frame the Thomas 
and the geodetic precession will vanish and we are left 
only with the Lense-Thirring precession. 

The spacetime to be considered in this work has the 
form, 



gn{dx l Y + g 2 2(dx 2 y + g 33 (dx 6 y + 2g 3i dx 6 dx 



+ 344 (dx 4 ) 2 . 

We will choose the frame associated to the tetrads 
19 15 
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y/-ga dx 4 ' 

where F = 333 — g 3i / ga (see also [ll|). The four velocity 
it M in this frame is given by u M = (0, 0, 0, (— 344)~ 1 ^ 2 )- 

To caculate the angular velocity of precession relative 
to the choosen tetrads we use the Fermi- Walker transport 
law, 



V u S = u(aS), a = V u u. 



Therefore 



dr 



V»(S.ej), 
S-(V u ej), 
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where V u is the total derivative operator, u the four- 
velocity and S the spin vector. Using the fact that 
ep,-B0 = r/^c and considering that the spinning particle is 
at rest relative the tetrads we have u = e 4 . Hence the 
condition u • S = gives us S% = 0. From Ve M = e v m v , 
we find that the Fermi-walker transport law can be cast 
as, 



dr 



(S-e S )(mVu), 
5 a (m a r e a ). 
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Now using the connection coefficients given in Appendix 
|A1 and omitting the subscript LT in Q.lt we find, 



— =flxS, 
dr 



where 



The symbols Tgji and are given by 
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We can write (fT6|) and (jXTJ) in a simpler way introducing 
the norm A and the twist oj^ of the time-like Killing vector 
field ^ = (0,0,0,1), 



A = 



(18) 
(19) 



From ©,(TTH1), OH) and ©-© we find that T m and T m 
can be written as, 



r _ 1 
m ~ 2X UJh 
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m ~ 2A 2 " 
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From the vacuum Einstein equations we know that — 
Hence the angular velocity can be cast as, 



V M w [1 
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The scalars A and u> are related to the Ernst potential, 
r, by the relation T = —A + iu> [T^l - The expression 
(|2ip is particularly interesting because it relates, in an 
exact manner, the angular velocity of precession to the 
two fundamental scalars of a stationary spacetime. We 
were unable to find this expression in the literature. 

To express the scalar ui in terms of Thorne moments 
using an harmonic coordinate system we use the defini- 



tion 
field 



and the Thorne metric for a weak gravitational 
]. This metric is given by, 



g^v = h MV + Tf MV , 7^ = - ^r]^ri al3 h a 0, (22) 
1=2 
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where 77^ is the Minkowski metric. Ma 1 and 5U, are the 
mass and angular momentum multipole moments respec- 
tively. A useful expression is, 



A = (-l)'(2l-l)!!*£. 



(26) 



where Nai — n ai ...n ai and n a — x a /r. 

Now from the previous equations we find that u> a can 
be written as 



From Eqs. 
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Now we shall write this solution in spherical coordinates. 
Since we are working in an axially symmetric spacetime 
our multipole moments M Al and S A , are multiples of 
z <Al> , the symmetric, trace-free outer product of the 
axis-vector with itself. The moments are completely de- 
termined by the numbers Mi and Si Q defined by 



Mi=M Al z A \ Si=S Al z A <. 
We can show that, 



(29) 
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To compute the quantity S Al N Al we use S Al N Al — 
Si Pi (cos 8), where P/(cos6>) are the usual Legendre poly- 
nomials. Hence, 



E 



41(21- l)l\SiPi(cost 
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To be consistent in our approximation we use (744 = — 1 
in OH). Then 



1. 



(32) 



Therefore in spherical coordinates, 

n = ± 2 -^^^[(l + l)Pi(cos8)f 



(1 + 1)! r l + 2 1 



sin 8P1 (cos 8) 8], (33) 



where ()' denotes derivative with respect to cos8. Equa- 
tion (|33p shows that, in this approximation, the contri- 
butions to are only from the angular momentum mo- 
ments. In Appendix [B] we show that the above formula 
is equivalent to formula (23) of [l4j. Relating i~2 to the 
Cartesian unit vectors x, y and z, we get 



^ 2i(2/-l)!! Si 

(/ + !)! ^ f( ^' 

where the components of f (I, 8) are 



(34) 



f x = (Z + l)P z (cos0)sin0 + P z (cos0)cos0sin0,(35) 
P = 0, (36) 
f z = (I + 1) Pi (cos 8) cos 8- P{ (cos 8) sin 2 8. (37) 

In deriving the previous expressions, for convenience, it 
was assumed that the orbit is on the plane y = 0. 



non-vanishing multipole moments this circle is slightly 
distorted. The radius of such orbit in the lowest order 
correction of the quadrupole moment is given in [TEI ] . In 
terms of Thorne quadrupole moment Q the orbit can be 
written as, 



r = r 1 - 



3Q cos26> 
8M rl 



(38) 



where M and Q are the mass and Thorne quadrupole 
moment of the source, respectively. From Eq. (I34p . ex- 
panding in powers of l/^o, we obtain 
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where the vectors gj(9), j — 1, 2, 3 are, 

4 16 
gl = 2f(l,0) g2 = -f(l,0)cos20, g 3 = yf(3, 



(40) 



The vectors gj have the property (g*) — (g^) = and 
(dj) = 1) where (g^) means to take the average over <?j (8) , 
which in this case is a simple integration over 8 and in 
view of the orbit symmetry, it is enough to consider half 
orbit from one pole to the other [0, it}. The average of ft 
reads, 



:n) = A 



(41) 



where 



A = 



S 

2rl 
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MSs 
QS 



The case C = 2/5 is special, since we have no correction 
although the spacetime is deformed. 

The Thorne mass moments can be accurately com- 
puted using de data published in [l6[, but the angular 
momentum moments are independent of the mass mo- 
ments. To determinate the constant C we need a model 
for an oblate Earth. Some authors [3 have already 
calculated the non-spherical contributions to fl. We shall 
use our approach to compare the precession for these 
models. Also, based on the Earth model of reference [l8[, 
we give a new estimative for ft. 



A. Teyssandier model 



III. n FOR THE GYROSCOPE EXPERIMENT 

The ideal orbit for the Gravity Probe B satellite is a 
circular one with altitude 642km. Since the Earth has 



Comparing the general form of Thorne metric with the 
one in [14|, we get 



Q = - 



2MR 2 J 2 



5, = 



4SR 2 K 2 



C = 0.97. (42) 
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where we have used the values J 2 = (1082.64 ±0.01) x 
10" 6 e K 2 = 0.874 x 10~ 3 given by the author. Therefore 
we will have 

(it) = A(l + \A2B)z. (43) 

This value is different from the one in because we 
use the angle 9 in p4|) instead of ip — ir/2 — 9 as in the 
quoted reference. Therefore on averaging over (|39|) we 
get a different, but equivalent result. 

B. Adler-Silbergeit model 

For Adler-Silbergeit model B fl7l |. we obtain 
2MB? J 2 



Q = — 



8ivMB 4 J 2 rAA . 
S 3 = — , S = Iu>. (44) 



From this expression we calculate the constant C = 
MSs/QS, and Eqs |@TJ gives us, 



in) =A 



(45) 



Confronting this result with Eq. (61) of [17J we see that 
our third term differs by a factor two (the numerical fac- 
tor in the quoted reference is 3/7). Using MB? 1 1 = 3.024 
given by the authors we obtain 



(ft) = A{1 + 1.59B)z. 
They get (ft) = A{1 + 0.30B)z. 



(46) 



Adler model 



For this model [18| we have 
2Ma 2 



Ss = - 4 4^, C = 0.72. 



9 ' 25 
Hence the averaged angular velocity is 

(Q) = A{1 + QMB)z 



(47) 



(48) 



where M n ,J n are the mass and current moments, re- 
spectively, and a, (3 and k are parameters. We find the 
following relation, 



J_ 

M 



(2M 2 



M Kerr 



(50) 



-J 2 /M is the Kerr quadrupole moment. 



where Mf err 

Since M 2 err has vanishing Newtonian limit, in our ap- 
proximation, we are only left with 

2M 2 J 



J3 = 



M 



(51) 



Now using the correspondence between Geroch-Hansen 
and Thorne moments [§] we get C — 4/15. Hence equa- 
tion (|4"T]) reads 

(n) = A(l - 0.3W)z. (52) 
that differs considerably from the previous models. 

IV. SUMMARY AND CONCLUSIONS 

By using an adequate tetrad basis we derive an exact 
expression for the Lense-Thirring precession in terms of 
the norm and the twist of the spacetime timelikc killing 
vector. We calculated the imaginary part of Ernst poten- 
tial using the weak field approximation and the Lcnsc- 
Thirring precession to any desired multipolar correction. 
We consider the case of the Gravity Probe B experiment 
expanding the precession angular velocity up to the or- 
der of 1/fQ and averaging over a half trajectory. The 
final expression is given in terms of Thorne multipole mo- 
ments. In our opinion, this general expression, may be 
interesting for the interpretation of exact solutions of Ein- 
stein field Equations. We have reproduced known results 
and given new estimatives for the non-spherical contri- 
butions to precession angular velocity. We conclude that 
the earth model plays a significant role in the multipolar 
contribution. The results produced by different Earth 
models are accounted for by our constant C, fact that 
makes our relation (|4"Tj) particularly usefull. 



D. Exact solution 

It is instructive to calculate the constant C for the case 
our geometry is described by an exact solution of the Ein- 
stein vacuum equations with arbitrary quadrupole mo- 
ment. We choose the version of this solution presented 
in [l!^ |. Equivalent solutions has been obtained by sev- 
eral authors using different methods [12]]. The first four 
nonzero Geroch-Hansen moments are, 

M = fc(l + a 2 )/(l-Q 2 ) 

J x = -2ak 2 (l + a 2 )/(l-a 2 ) 2 

M 2 = -fc 3 [/3 + 4a 2 (l + a 2 )(l - a 2 )" 3 ] (49) 

J 3 = 4a/fc 4 [/3 + 2a 2 (l + a 2 )(l-a 2 r 3 ]/(l-a 2 ) 



Acknowledgments 

The work of MZ was supported by Capes. PSL thanks 
CNPQ and FAPESP for partial financial support. 

APPENDIX A: CONNECTION COEFFICIENTS 

For easy reference we give in this appendix the connec- 
tion coefficients relative to the tetrads ©-(HJ). The dual 
base is: 

m = ^gudx 1 , m 2 = ^g 22 dx 2 , (Al) 
m 3 = y/Fdx 3 , m 1 = ; 934 - ^g^ldx A .{A2) 

V~94i 
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From the first Cartan structure equations [l2j |. 

dm'' + mt x ra v = 0, 
we find the connection coefficients, 



where 



and 
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APPENDIX B: EQUIVALENT EXPRESSIONS 

In this appendix we show that our Eq. (|33[) is equiva- 
lent to Eq . (23) of Ref. [14} . Using the relation, 



sin f 



r(cOS#f — z), 



(Bl) 



we can write equation f|33[) in terms of the vectors r and z, 



n ^ 2i(2Z-l)!! Si „, ... „. 



+ cos 6»P ; (cos 9)f-P l (cos 0)z]. (B2) 

From the relation, 

P ; ' +1 (cos6») = (/ + l)P;(cos0) + cos 6>P,' (cos 0), (B3) 
we can cast (|B2[) in the form, 



(B4) 



or as, 



ft = — (3cos0f — z+ 
^(^^'^[^(cos^-P/Ccos^). 



^ (/ + 2)Z! r 



(B5) 



Identifying, 

= {l + 2)\SKiR l 
1+1 2(2/ + 1)!! ' 

we get Eq. (23) of Teyssandier paper. 
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